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Abstract 

We study the processes e'^e' — > (J/'i/'77)7 and e^e^ {J / ■iIjtt~ Tr^)'y where the hard photon 
radiated from an initial e^e~ collision with center-of-mass (CM) energy near 10.58 GeV is detected. 
In the final state J/ip^^ we consider J/^'vr'^, J/ipr], Xcil-, and Xc2l candidates. The invariant mass 
of the hadronic final state defines the effective e^e~ CM energy in each event, so these data can 
be compared with direct e^e~ measurements. We report 90% CL upper limits for the integrated 
cross section times branching fractions of the J/'i/'77 channels in the y(4260) mass region. 
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1 INTRODUCTION 



In 2005, BABAR reported the first observation of a broad structure in the J /Tp'n~ii~^ mass spectrum 
produced via e'^e~ annihilation and detected in initial state radiation (ISR) events: e'^e~ — > 
{J / 'tpn~ TT~^)^ . The mass was reported to be approximately 4.26 GeV/c^ and the width ~ 90 MeV/c^ 
if interpreted as a single resonance [1]. Subsequently, CLEO-c reported cross section measurements 
for J/'ip'K~Tr~^ and J/ipir^ir^ production at this mass, confirming the large production rate [2]. More 
recently, they have reported ISR measurements of the J / ipn" n'^ final state [3]. As discussed in refs. 
[1], [2], [3], and references therein, the nature of the y(4260) is not clear, so observing additional 
decay modes and measuring their relative rates is important for advancing our understanding of 
this state and its production rate. 

In this paper we report measurements of e^e~ — > J/'077 cross sections using ISR events where 
the hard photon radiated from an initial e~^e~ collision is detected directly. Requiring that the ISR 
photon is detected reduces backgrounds significantly and increases the signal-to-background ratio, 
which we confirm with a parallel study of the benchmark signal e"^e~ — > J / tpn" tt'^ . In the final 
state J/ip'yy we consider J/ipir^, J/'ipV^ Xcil, and Xc27 candidates. If the y(4260) is a charmonium 
state, its decays might be similar to those of the ^'(^S) which has relatively large branching fractions 
to J/ipr], Xcil, and Xc2l [4]. 

The ISR cross section for a particular hadronic final state / (excluding the radiated photon) is 
related to the corresponding e+e" cross section crf{s) by: 

^^I^ = W{s,x)-aM^-^))^ (1) 

where x = 2E^l^fs\ is the energy of the ISR photon in the nominal e+e" center-of-mass (CM) 
frame; y/s is the nominal e'^e~ CM energy; and \/ ,5(1 — x) is the effective CM energy at which the 
final state / is produced. The invariant mass of the hadronic final state defines the effective e"'"e~ 
CM energy. The function W{s, x) is calculated with better than 1% accuracy [5] and describes the 
probability density function for ISR photon emission, which occurs at all angles. For the present 
study we require that the hard ISR photon be detected in the electromagnetic calorimeter of the 
BABAR detector. The effective ISR luminosities at CM energies corresponding to the masses of the 
J/ifj, the -0(25), and the y(4260) are 66.3 MeV^^ nb"\ 84.3 MeV"^ nb"\ and 105 MeV^^ nb'^ 



2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The data set used here is the same data set studied in the original y(4260) analysis [1]; the 
integrated luminosity at and near 10.58 GeV in the e~^e~ CM is approximately 230 fb~^. 

Charged-particle momenta are measured in a tracking system consisting of a five-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer central drift chamber (DCH), both situated in a 
1.5-T axial magnetic field. An internally reflecting ring-imaging Cherenkov detector (DIRC) with 
bar radiators made of fused silica provides charged-particle identification. A Csl electromagnetic 
calorimeter (EMC) is used to detect and identify photons and electrons. Muons are identified in 
the instrumented magnetic flux return system (IFR). 

Electron candidates are identified by the ratio of the shower energy deposited in the EMC to 
the momentum, the shower shape, the specific ionization in the DCH, and the Cherenkov angle 
measured by the DIRC. Muons are identified by the depth of penetration into the IFR, the IFR 
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cluster geometry, and the energy deposited in the EMC. Pion candidates are selected based on a 
likelihood calculated from the specific ionization in the DCH and SVT, and the Cherenkov angle 
measured in the DIRC. Photon candidates are identified with clusters in the EMC that have a 
shape consistent with an electromagnetic shower but without an associated charged track. A 
detailed description of the BABAR detector is available elsewhere [6]. 

We use an admixture of simulation packages for the final states studied. These packages use 
adaptations of algorithms for generating hadronic final states via initial state radiation [7], for 
modeling multiple soft-photons from initial-state electrons and positrons via a structure-function 
technique [8, 9], and for modeling photons from final-state particles [10]. We pass the generated 
events through a detector simulation based on GEANT4 [11], and reconstruct them in the same 
way as we do the data. 

3 EVENT SELECTION CRITERIA 

The event selection criteria used in this study, and described below, were chosen to produce very 
low backgrounds and to eliminate events with relatively poor mass resolution. These goals are 
especially important for the high-mass J/ tp'j'j final state where the anticipated signal level is small. 
We therefore require that the ISR photon be detected directly for this analysis and do not use the 
much larger sample of events where the ISR photon is produced outside the EMC acceptance. To 
avoid bias, possible signals in the mass region above 3.85 GeV/c^, were not studied until all of the 
selection criteria had been determined. 

All charged tracks used in this analysis are required to have laboratory polar angle in the 
range 0.45 to 2.40 radians and are required to have at least one hit in the SVT and at least 
20 hits in the DCH. Such tracks are reconstructed relatively well and are detected in a fiducial 
region where the particle identification is understood especially well. Candidate J/V' mesons for 
the J/'i/'77 analysis are identified via their decays to ^'^ ^~ . Those for the J/ipTr~Tr~^ analysis are 
identified via their decays to and e+e~. We require that both tracks in a dimuon candidate 

be identified as muons, with at least one satisfying a set of criteria which may be referred to as 
loose and the other satisfying a set of criteria which may be referred to as very loose. Using a 
sample of ip{2S) — > J/ipir~7r~^ decays produced in ISR events, we measure the muon identification 
efficiency for pairs to be f« 80%. Because QED production of four-lepton events is a major source of 
background for the J/V'vr~7r+ sample, we exclude all candidates where either of the pion candidates 
is identified as possibly a muon or an electron. To improve the signal-to-background ratio, and 
because the J/V' — > e~^e~ signal has a long radiative tail and thus spans a greater mass range 
than does the J /ip ^ signal, we require that both charged tracks be positively identified as 

electrons in selecting e^e~ pairs. We do not measure the electron pair efficiency per se, but measure 
the J/^p e+e" efficiency relative to the J/^p /^"""M" efficiency using ISR ^(2S) J/ip-K~Tr^ 
events. 

ISR events are characterized by small missing mass, corresponding to that of the ISR pho- 
ton, recoiling from the exclusive hadronic final state of interest. Experimentally, the missing mass 
squared, measured using the beam momentum and the momenta of the final state hadrons (includ- 
ing the non-ISR photons), is not exactly zero. Most importantly, higher order QED contributions 
to ISR production of the final state produce a very long tail at positive missing mass. In addition, 
the finite spread of the beam energy/momentum and the the detector resolution further distort the 
observed values. However, the signal is predominantly produced with small recoil mass, and the 
signal-to-background is greatest when the observed recoil mass squared is close to zero. We exploit 
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this feature of the process both to select a sample of candidates with lower background rate and 
to improve the resolution of events selected. We require the recoil mass squared, determined from 
the final state particles' measured momenta and the beam particles' nominal momenta, to lie in 
the range -2.0 to 4.0 GeV^/c^. We also require that the missing energy, calculated in the lab, fall 
in the range -0.4 GeV to 1.5 GeV and that no more than 400 MeV of additional neutral energy be 
detected in the EMC. We then re- fit the hadronic final state particles' momenta using the nominal 
beam particles' momenta, the relevant track parameter error matrices, and the constraint that the 
recoil mass is zero. We call this the IC fit and require its to be less than 20. The IC fit predicts 
a direction for the ISR photon. We require that the direction of the observed ISR photon candidate 
coincide with this predicted direction within 15 mr. In addition to reducing backgrounds, this last 
requirement, in the absense of any others, rejects approximately 25% of all signal candidates as a 
result of the higher order QED processes that contribute to the ISR cross section. In such events, 
the IC fit must also be high, so the pointing and criteria are strongly correlated. 

We select as J /ip candidates dimuon pairs whose mass lies within 20 MeV/c^ of the nominal 
J/ip mass [4]. The IC fit improves the dimuon resolution significantly, but as our charged track 
resolution is better than our photon resolution, this effect is less pronounced for the J/ip^^ sample 
(where the J /il) mass resolution is ~ 12 MeV/c^) than for the J sample (where the J/V' mass 
resolution is ~ 7 MeV/c^). However, as the signal-to-background (in the sample of J/V' candidates) 
is worse in the J/'i/'77 sample, we use the same dimuon mass criteria for both J/iI:ii~'k~^ and J 
final states. For the J/ipii~T:~^ channel we also select as J/ip candidates electron pairs whose mass 
lies in the range 3.050 - 3.120 GeV/c^. This asymmetric window around the nominal J/il) mass 
is used because the signal has a long radiative tail. The specific limits of the window were chosen 
following studies of ISR production of 1/^(28) — > J/ip'K~TT'^ with J/ip — > e"^e~. We do not use 
J/'0 ^ e^e~ candidates in this study of J j^)^^ final states as the ISR data set used in this analysis 
was defined before the discovery of the y(4260), and excluded e+e~777 candidates intentionally. 

For the 7/^^77 candidates accepted by the selection criteria described above, a two-constraint 
(2C) fit is done where the invariant mass of the lepton pair is additionally constrained to be 
that of the J/V'- We then explicitly require the presence of intermediate states so that J/'i/'77 
is a manifestation of J/ijjTi^ ^ Xcil, or Xc27- We require that each of the non-ISR photons 

has measured energy in the laboratory of at least 100 MeV and that its laboratory polar angle 
lies in the range 0.35 to 2.40 radians, where we understand the EMC performance well. Using 
momenta from the 2C fit, we require that the invariant mass of a photon pair lies in the range 
0.120 - 0.150 MeV/c^ to be considered as a vr candidate and in the range 0.52 - 0.58 MeV/c 
to be considered as an r/ candidate. Similarly, we require that a J/'i/'7 invariant mass lies in the 
range 3.4906 - 3.5309 GeV/c^ (3.5363 - 3.5763 GeV/c^) to be considered as a Xci (Xc2) candidate. 
Events with such candidates are then re-fit with a third constraint requiring that the invariant 
mass of the corresponding combination be that of the hypothesis. When an event satisfies the 
nominal requirements for more than one hypothesis, the hypothesis with the best 3C x^ is selected 
as that to be used in the final analysis. Given our 77 and J/'i/'7 resolutions, these windows accept 
about 90% of the candidates otherwise accepted. The J/V'77 resolutions depend upon the mass 
of the final state, in particular how far above threshold it is. For the decay y(4260) J/il'iJ^ 
(simulated) resolution is ~ 10 MeV/c^. For all of the y(4260) decay modes considered, the mass 
resolution is negligible relative to its ~ 90 MeV/c^ width. 
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Figure 1: The upper plot shows the J/il^r] 
invariant mass for J/ij}^^ candidates where 
J/ipr] is the best hypothesis. The lower plot 
shows the iI)tt^ invariant mass for J jil)^^ 
candidates where J /ij)T:^ is the best hypoth- 



Figure 2: The upper plot shows the Xci 7 
invariant mass for J/ip^^ candidates where 
Xci 7 is the best hypothesis. The lower plot 
shows the Xc2 7 invariant mass for J/ 1/;77 
candidates where Xc2 7 is the best hypothe- 
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4 MEASUREMENTS AT THE V^(2S) 

We initially look at signals in the '0(2S) mass region but not above 3.85 GeV/c^. We also estimate 
background levels by looking at candidates with dimuon invariant mass in J/ip sidebands. The 
■0(28) mass region J/ijj^^ invariant mass distributions for candidates selected as J/ipr] or J/ip-K^ 
are shown in Fig. 1. The corresponding distributions for those selected as Xci7 or Xc27 are shown 
in Fig. 2. Note that the mass ranges of the samples differ. In particular, the J/tpn^ threshold is 
much lower, and therefore its mass resolution is much worse. 

The signal levels are low in all four channels, and the background levels much lower, so we 
estimate the signal in each channel as the number of entries within 15 MeV/c^ of the nominal V'(2S) 
mass. This gives 46, 4, 84, and 34 events for the J/il^rj, J/ipT^^, Xcil, and Xc27 channels. Monte 
Carlo simulations indicate that the efficiencies for these channels are roughly equal (not including 
the branching fractions of the intermediate states). The probability that an event generated as Xci7 
will be reconstructed as Xc27 or as J /ipr] is of order 5%. The probabilities that an event generated 
as J/il^rj or J/ipir^ will be reconstructed as Xci7 or Xc27 are of order 5% each. Altogether, the 
number of ^p{2S) J/i^H signal events we observe is generally consistent with, but somewhat 
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higher than, the number projected from Monte Carlo studies using branching fractions (or products 
of branching fractions) from CLEO [12], which are generally higher than those reported by the PDG 
[4] in 2004. 

For both samples of J I'^tx~tx^ candidates accepted by the selection criteria described above, a 
two-constraint (2C) fit is done where the invariant mass of the lepton pair is additionally constrained 
to be that of the J/V'- We observe clean "0(28) signals with 785 candidates in the J z""^/^" 
sample and 434 in the J/ip ^ e^e~ sample. In both channels the width of the peak, fit as a 
Gaussian, is « 3 MeV/c^ and the background rate is at the percent level. From the relative 
numbers of ■0(28) signal events in the two J/V' decay modes, we find the combined acceptance 
times reconstruction efficiency to be 1.55 times that for the J/ijj ^ A^"*"/^" mode alone. We assume 
the same ratio holds at higher J/iIjtt~tt^ mass. 

The integrated ISR luminosity, requiring the ISR photon to be produced in the laboratory with 
a polar angle in the range 0.35 to 2.40 radians, where we understand the EMC performance well, 
is 84.3 MeV"^ nb~^ at the mass of the '(/^(2S), and we estimate its fractional uncertainty to be 3%. 
For a narrow resonance, one may write the integrated cross section 



^ = j o'(Ecm)'^Ecm 



where a is a resonant cross-section and is assumed to have a Breit-Wigner line shape. It is directly 
related to the partial decay rate to electron pairs. If the resonance has spin 1, then 

i:=^ree. (3) 

Using PDG values for Fee, rn, and relevant branching fractions [4], we predict 784 "0(28) 
J/ip'K~'K'^] J/ij} ^ fJ'~^fJ^~ events with an 8% fractional error and we observe 785. In a similar 
exercise, we predict ^ 17,900 ISR J/^p — > iJ,~^fJ-~ events with a 4.6% fractional uncertainty and we 
observe 17, 312 it 271 events, about 97% of the predicted number. From these studies we conclude 
that we understand the product of ISR luminosity and J/'07r~7r~'^ reconstruction efficiency with no 
worse than 10% precision. 



5 MEASUREMENTS AT HIGHER MASS 

We estimate high-mass J/'077 background rates due to non-J/ip dimuon pairs in the J/^l' mass 
range using fi~^fJ.~ background windows 75 to 175 MeV/c^ above and below the nominal J/tp mass. 
For dimuon pairs in these windows we look for rj and vr'' candidates as we do for real data (except 
for performing the 2C fit). To create a sample of xc candidates, we calculate the /i~^"^7 invariant 
masses, subtract the /^"/i''' mass, and add the J ftp mass to compare to the Xci and Xc2 masses. 
In all four cases, we calculate an effective "J/'077" mass as the /i~//~''77 mass minus the 
mass plus the J/tp mass. These background windows enclose a total of 200 MeV/c^ while the J/ip 
candidates lie within a 40 MeV/c^ wide window, so the background is estimated by dividing the 
number of candidates observed by 5. This predicts between 2 and 3 background events in the J/tprj 
sample, predominantly found below 5 GeV/c^, and about half this amount in each of the other 
channels. 

The high mass J/tpr] and J/tpir^ data are shown in Fig. 3. The high mass Xci 7 and Xc2 7 
data are shown in Fig. 4. We see no events in the y(4260) region in the J/ipr], J/ipi:^, or Xc2 7 
distributions, and about the total number of entries estimated as non-J/'0 fi~^fi~-y^. Two entries 
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Figure 3: These are the J/il^rj (top) and Figure 4: These are the Xcil (top) and 
J/V'tt" (bottom) invariant mass distribu- Xc2 7 (bottom) invariant mass distributions 
tions for candidates satisfying all final anal- for candidates satisfying all final analysis cri- 
ysis criteria. teria. 

are observed in the y(4260) region of the Xci 7 distribution, and six entries altogether, about two 
or three times the number estimated as non-J/'0 /i~/i^77. This leads to the question of other 
potential sources of background. 

Excluding J/ijj^^ candidates whose invariant mass lies within 25 MeV/c^ of the ■0(28), we 
examined the remaining 77 invariant mass distribution and separated the pairs into those from each 
of the four candidate categories and those in which the parent did not satisfy the criteria for any of 
the intermediate states. We use this last set of events to estimate how many additional background 
events might be observed in the Xci7 sample in the y(4260) mass region. The statistics of the 
samples are low, so a precise estimate is not possible. A reasonable estimate is that the expected 
background rate is between 0.5 and 1.0 additional events. Thus, we cannot exclude the possibility 
that both events observed are background. Neither can we exclude the possibility that they are 
signal. We therefore calculate the 90% confidence level upper limit assuming no background. 

We calculate the 90% confidence level upper limit for the integrated cross section times branch- 
ing fraction (or product of branching fractions) for each of the four channels considered using 2.3 
events as the Poisson upper limit for channels where no entries are observed in the y(4260) region 
and 5.3 events in the case of the xci-f channel where 2 entries are observed. Similarly, we use 5.3 
events for calculating the upper limit for inclusive J/'077 production of these states. The Monte 
Carlo efficiencies for the J/tpr], J/ipir^, and xc7 channels are 11.2%, 10.9%, and 11.9% respectively. 
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The efficiency for inclusive production of these states, 12.2%, is somewhat higher than the 

efficiencies of the exclusive final states because, in contrast to the those efficiencies, it counts as 
reconstructed signal also those events generated as one exclusive final state but reconstructed as an- 
other. To allow for systematic uncertainties in luminosity times reconstruction efficiency, we round 
off each of the precisely calculated values, increasing them by between 10% and 20%. These results 
are shown in Table 5. The ratio of branching fractions, with respect to y(4260) — > J/i/;7r~7r"'", 
is calculated using the cross section times branching fraction result from our original study [1]: 
H X e(y(4260) ^ J/?/'7r-7r+)) = (7.0 ± 1.3 ± 1.0) MeV nb. 

Table 1: This table compares the integrated cross section times branching fraction results for the 
y(4260) — > J/'i/'77 decay modes studied. Limits relative to y(4260) — > J/iI)t^~tt~^ are calculated 
using the cross section measured in our original study [1]. 



channel {X) 


E X B(X) 


B{X)/B{J/^p^^-^^+) 


y(4260) ^ J/V'r/ 


< 10 MeV nb 


< 1.4 


y(4260) ^ J/V'vr" 


< 4 MeV nb 


< 0.6 


y(4260) ^ Xcil ^ J 1^11 


< 25 MeV nb 


< 3.6 


y(4260) ^ Xc27 ^ ^/V'77 


< 18 MeV nb 


< 2.6 


y(4260) ^ J/ipll 


< 8 MeV nb 


< 1.2 



As it is our benchmark mode, we also study J/tpn'TT'^ at high mass using selection criteria 
very similar to those used for J/ip^^. We first studied expected background rates. We defined 
wide dilepton mass windows above and below the J /ip mass region in both channels and used the 
^"'"£~7r^7r~ candidates in those windows as if they were J/ipTT~TT~^ candidates to estimate background 
rates. Scaling the widths of the background windows to those used for J/ip candidates, we estimate 
8 background events in the 2.5 GeV/c^ range from 3.75 GeV/c^ to 6.25 GeV/c^, distributed more or 
less uniformly over that range and divided roughly equally between J /ip ^ and J /ip ^ e^e~ 

candidates. 

The high mass J/'i/'vr~7r+ sample for this study is shown in Fig. 5. For comparison, the plot 
from our original y(4260) paper [1] is shown in Fig. 6. As anticipated, the signal level is much lower 
in this study and the signal-to-background ratio is much higher, primarily because here we require 
that the ISR photon be detected directly and not produced outside the EMC acceptance. The 
non- J ftp /i~/i"''7r~7r"'" background rate in this plot, estimated from the J/tp sidebands, is « 0.064 
events per 20 MeV/c^ bin. It is insufficient to explain the 10 entries observed in the mass range 
5.0 GeV/c^ to 6.25 GeV/c^, where it predicts only 4. In the window ±250 MeV/c^ around 4.260 
GeV/cP we observe 24 events where our background prediction is 1.6. The data are consistent with 
a structure 90 MeV/c^ in width, but the statistics are too low to draw an independent conclusion. 
In the original discovery of the y(4260) [1], the resonant cross section times branching ratio was 
measured to be (7.0 it 1.3 it 1.0) MeV nb. This predicts a central value for the resonant signal level 
in the current analysis to be ~ 10.5 it 2.5 events. The higher event rates observed in the y(4260) 
mass region and at higher mass suggests the possibility of continuum J/tpTT~7r~^ production. 
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Figure 5: This is the J/^^-tt+ invariant ^(^260) paper [1]. 

mass distribution for candidates satisfying 

all final analysis criteria. 

6 SUMMARY 

We observe ISR production of ^(2S) decaying into J/il^rj, Xci7, and Xc27 with summed branching 
fraction generally consistent with, but somewhat greater than, that reported by CLEO. We see 
no events in the y(4260) region in the J/ipi], J/ipir^, or Xc2 7 distributions. We observe two 
entries in the y(4260) region of the Xci 7 distribution, which is consistent with being either signal 
or background. We set 90% confidence level upper limits on the integrated cross section times 
branching fraction (or product of branching fractions) for each of these channels. 

We observe ISR production of J /ipT:~iT~^ events near 4260 MeV/c^ consistent with the mass 
and width originally reported for the y(4260) [1]. With the improved signal-to-background ratio 
of the analysis reported here, we observe an excess of events, both in the y(4260) mass region and 
at higher mass. These observations suggest continuum production of J/ipn~7r~^. 
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